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We present an analysis ol the mean-field kinetics of Brownian coagulation of droplets and polymers 
driven by input of monomers which aims to characterize the long time behavior of the cluster size 
distribution as a function of the inverse fractal dimension, a, of the aggregates. We find that 
two types of long time behavior are possible. For < a < | the size distribution reaches a 
stationary state with a power law distribution of cluster sizes having exponent §. The amplitude 
of this stationary state is determined exactly as a function of a. For | < a < 1, the cluster size 
distribution never reaches a stationary state. Instead a bimodal distribution is formed in which a 
narrow population of small clusters near the monomer scale is separated by a gap (where the cluster 
size distribution is effectively zero) from a population of large clusters which continue to grow for all 
time by absorbing small clusters. The marginal case, a — |, is difficult to analyze definitively, but 
we argue that the cluster size distribution becomes stationary and there is a logarithmic correction 
to the algebraic tail. 



I. INTRODUCTION 

Aggregation underlies numerous phenomena from milk 
curdling and blood coagulation to planet and star forma- 
tion. A review of the earlier work on aggregation has 
been given by Chandrasckhar [l|; the following work is 
summarized in a number of books and reviews [2-8]. This 
research field has been initiated by a pioneering work of 
Smoluchowski [9| who established the framework for the 
analysis of aggregation. Smoluchowski considered well- 
mixed diluted systems and argued that the governing 
equation for the concentration c/. (t) of clusters of "mass" 
k (that is, clusters composed of k monomers, where the 
monomers are the clusters of the minimal mass) reads 



dc k 
dt 



1 
2 



E 



K, 



T.K, 



i.jCiCj Cfc 2 , 

3>1 



k,j c j 



(1) 



Here Kj is the rate at which clusters of mass i and j 
merge. Only binary collisions are taken into account since 
the system is assumed to be diluted. Mathematically, the 
reaction rates Kij form an infinite symmetric matrix, 
Ki.j = Kj,i\ the pre- factor ^ in front of the gain terms 
on the right-hand side of (flj is needed to avoid double 
counting. 

In addition to devising the mathematical framework 
underlying aggregation processes, Smoluchowski also 
computed the reaction rate Kij in the most important 
case of Brownian coagulation. Namely, assuming that 
aggregates are spherical, Smoluchowski showed that 



K itj = 47r(A + DjXRi + Rj) 



(2) 
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where D^ and Rk are the diffusion coefficient and the 
radius of the cluster of mass k. Invoking the Stokes- 
Einstein relation between the diffusivity of a spherical 
(three-dimensional) object and its radius yields D ~ R^ 1 
and hence Eq. © gives Kij ~ (R^ 1 + Rj 1 )(R i + Rj). 
Using additionally volume conservation, Rk ~ fc 1 ' 3 , one 
arrives (up to a multiplicative factor that can be absorbed 
into the time variable) at the following expression for the 
Brownian reaction rate 

^=Q) 1/3 + (f) 1/3 + 2 (3) 

Rate equations (JIJ with the Brownian reaction rates 
§3§ have never been solved. It is not shocking, of course, 
as mathematically they form an infinite set of coupled 
non-linear differential equations (JT|) with varying coeffi- 
cients (|3]). For simpler reaction rates, Eqs. (fT]) can admit 
analytical solutions. For instance, Smoluchowski already 
noticed Q that the Brownian reaction rate is homoge- 
neous, K\i,\j = Kij, This feature suggests to consider 
the simplest such reaction rate, Kj = const. In this sit- 
uation the rate equations ((T|) are indeed solvable 0, Q . 
A very few other solvable cases have been found in later 
work (see 0, |g-Q for review) , yet for all these years there 
has been no progress for the Brownian reaction rate. This 
classical case appears as analytically intractable problem 
as it originally looked. 

In this work, we consider aggregation with input. 
Such driven aggregating systems often approach a non- 
equilibrium steady state [l(j-[l2[ and the steady states 
tend to be more tractable. The details of input play 
rather limited role, e.g. they do not affect the emerging 
large mass behavior, if only clusters of small mass are 
injected. It is customary to assume that only monomers 
are injected. The governing equations then read 

" n l ^ - ^' "Ax (4) 



dt 



jCjCq 



2 / j &i,3^' f i 
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where we set the strength of the monomer flux to unity. 
(The flux strength can always be absorbed into the time 
variable.) 

In the long time limit, the mass distribution can be- 
come stationary. In this situation one must solve 



i+j=k 



i('i( ' i 






j^j 



5fc,l 



(5) 



The stationary solutions, when they exist, often exhibit 
power-law behavior. These are non-equilibrium station- 
ary states which arise due to the flux of mass from small 
mass scales to large mass scales supplemented by input 
at the small mass scale. These states are analogous to 
the flux-dominated stationary states that occur in fluid 
turbulence [13J . passive scalar advection [14J], wave tur- 
bulence pjl , granular gases [16| , and other driven aggre- 
gation systems fl7H2fj |. 

In this paper we present a detailed analysis of the long- 
time behavior of the cluster size distribution for range of 
physically relevant models of Brownian coagulation with 
input of monomers which are parameterized by a single 
scaling exponent, a, which is related to the fractal di- 
mension of the aggregates. We start by describing the 
parameter space of relevant models in Sec. |TT] and out- 
line where problems of Brownian coagulation fit in this 
space. Next in Sec. UTT] we study the cases < a < i 
which includes Smoluchowski's original problem of Brow- 
nian coagulation of spherical droplets. We show that, in 
this regime, the size distribution becomes stationary for 
large times. For the case of spherical droplets, we show 
that the tail of the size distribution has the form 



Cfc -> 



Ck^' 2 . 



C 



(6) 



in the k — > oo limit, the first exact result for this classical 
problem to best of our knowledge. In Sec. IIVI we study 
the cases 0| < a < 1 where we find that the size distri- 
bution does not become stationary and characterize the 
dynamics for such kernels which include the physically 
interesting cases of stiff polymers (a — I) and polymers 
in an ideal solvent (a ~ |). Finally in Sec. El we provide 
some insights into the marginal case of an ideal polymer, 
a = | although a complete analysis eludes us. We close 
with a short discussion. 



II. MODELS OF BROWNIAN COAGULATION 

Before beginning our study Brownian coagulation of 
polymers, let us first briefly describe the wider parameter 
space of models in which Brownian coagulation sits. This 
will help to place our subsequent work in context. 



A. A two-parameter class of reaction rates 

The two-parameter class of homogeneous reaction 
rates given by 



^,, = *v+*y 



(7) 



has been extensivelystudied (mostly without input) in 
the literature. See [7| for a review. The reason is that 
members of this class, for different values of the param- 
eters a and b, reproduce various rates (or their asymp- 
totics) that appear in applications. 

We shall only discuss models for which the exponents 
a and b satisfy 



b< 1. 



(8) 



The reason is that a system with kernel ([7]) undergoes 
a gelation transition when a + b > 1. Gelation corre- 
sponds to a loss of mass conservation in Eq. ([lj after a 



finite time [21[. In the absence of a source of monomers, 
a gelling system reaches a steady state, but the steady 
state is trivial, namely the entire system is comprised of 
a single cluster which, in the limit of an infinite system, 
has infinite size. (In the presence of a source, a station- 
ary state can be reached for such gelling systems [19| but 
they contain only finite total mass which considerably 
complicates the question of how the stationary state is 
approached dynamically. In particular, the scaling anal- 
ysis outlined in the appendix fails since the first moment 
of the scaling function diverges.) 

In addition, following the majority of work in this field, 
we shall further restrict the possible values of the param- 
eters a and b to the range 



a< 1, 



b< 1. 



(9) 



The mathematical reason is that when a > 1, and/or 
b > I , the system exhibits a highly singular phenomenon 
known as instantaneous gelation (see [8]). This essen- 
tially means that the solution of the infinite set of rate 
equations is ill-defined. One can make sense of a regu- 
larized version of Eq. CD in this regime but the cut-off 
plays a crucial role [22|]. In any case, for systems where 
coagulation requires physical contact between clusters, 
the reaction rate cannot grow faster than the masses of 
the clusters. The requirements © are then natural in 
such cases. 

It is instructive to recall the derivation [12| of the sta- 
tionary state for the rates ([7]), modulo the restrictions 
on the values of a and b discussed above. We shall em- 
ploy the generating function approach [2J|. We denote 
generating functions by 



e Q (z) 



k>l 



CfcZ 



and the corresponding moments by 



(10) 



(11) 



fe>l 



To find a stationary solution to ((4]), ([7]) we multiply Q 
by z k and sum over k to obtain 



e a (z)e b ( z ) - e a (z)M b - e b ( z )M a + z = o. 



(12) 



Let us assume that the cluster size distribution decays al- 
gebraically for small cluster sizes. Various exact solutions 
for particular kernels and the scaling argument outlined 
in Appendix [Al suggest that algebraic decay is plausible: 



c& 



C 



when k ^$> 1 



(13) 



To determine the decay exponent r and the amplitude 
C let us recast the conjectural asymptotic behavior of 
the mass distribution into a singular behavior of the gen- 
erating function. Indeed, the algebraic large k behavior 
(|13p of Ck is equivalent to the following singular behavior 
of the generating functions G a (z), Q b (z) in the z t 1 limit: 



e a = M a + cr(i - r + a) (i - zy-*- 1 

e 5 = M b + CT(1 - r + 6) (1 - zf- b - 1 H 



(14a) 
(14b) 



Substituting expansions (|14al) - (|14b|) into equation (fT2|) 
and matching the constant terms we obtain M a M b = 1. 
Matching then the leading (1 — z) behaviors we recover 
known results [12| for the decay exponent 



3 + a + b 



(15) 



An expression for the amplitude is rather cumbersome, 
but using standard identities for the Euler's Gamma 
function it simplifies to 



C 



[1 - (a - 6) 2 ] cos[7r(a - b)/2] 



47T 



(16) 



Using Eqs. (rT5"|) - (fT6"]) we can determine when the sys- 
tem can reach a steady state. First, we must assure that 
the series constituting M a and M b are convergent. This 
implies that r > 1 + a and r > 1 + 6, which in conjunction 
with ([15]) lead to 



\a-b\ < 1 



(17) 



This requirement simultaneously guarantees that the ex- 
pression (ITB)) for the amplitude is acceptable. The valid- 
ity of (JTTJ) and © actually ensures the validity of ©• 
Geometrically the region in the (a, b) plane where (fT7|) 
and ([5]) are satisfied, and thus a steady state is reached, 
is the semi-infinite strip denoted in Fig. [T]by region A. 

It is possible that the steady state is never reached. 
The analysis above shows that this can happen only out- 
side the semi-infinite strip denoted by region C in Fig. [T] 
For instance, a sub-monolaycr epitaxial growth is de- 
scribed by the generalized sum-kernel reaction rates 



Ki 



+ r 



(18) 



These rates are included in the class of rates ([7]) if one sets 
a = —jjL and 6 = 0. Two physically interesting cases are 




FIG. 1. Overview of the parameter space of the rate equa- 
tions Q with reaction rates (0. Kernels in region A admit 
a steady-state solution. Kernels in region B admit a steady 
state solution which is preceded by a gelation transition. Ker- 
nels in region C do not admit a steady state solution. Kernels 
outside of these regions undergo instantaneous gelation mean- 
ing that the rate equations are ill-defined in the absence of a 
cut-off. The solid line b = —a correspond to the family of 
Brownian coagulation kernels studied in detail in this paper. 



/i = 1 (this occurs when ad-atoms undergo "terrace" hop- 
ping) or /i — | (for the "periphery" hopping mechanism) 
|24| . In both of these cases, and generally when [i > 1, 
the system evolves ad infinitum [18], e.g. the cluster den- 
sity grows as (lni)^/ 2 for /i > 1. In the borderline case of 
ji = 1 , there appears an additional nested logarithm [18| , 
viz. the cluster density increases as ^/(lnt)/[ln(lnt)]. 
These behaviors are difficult to confirm by solving nu- 
merically mean-field rate equations, let alone to observe 
in simulations of a sub- monolayer epitaxial growth [251 ] . 



B. Brownian coagulation of polymers 

Let us now restrict our attention to the main topic 
of this paper, the rates for Brownian coagulation. The 
generalization of the Smoluchowski formula, @, from 
three to d dimensions is Kij ~ (Di + Dj)(Ri + Rj) d ~ 2 . 
In the physically interesting two- and three-dimensional 
settings, the range of relevant models is rather limited. 

If aggregation occurs on a two-dimensional substrate 
[26J . then K hJ - (D t + Dj). Now the diffusion coeffi- 
cient often varies algebraically with the size of the clus- 
ter, Dk on k~^ with /i being the mobility exponent, and 
this leads to the kernel (J18I) ; the corresponding behaviors 
have been investigated in Refs. |18j . 

In three dimensions for spherical aggregates we arrive 
at ([3]) . For fractal aggregates one would get a generalized 



Brownian kernel 



Ki 



0.25 



(19) 



where a = 1/Dj is the inverse fractal dimension of clus- 
ters. Most of the rest of this paper is devoted to the study 
of the solutions of Eq. ((TJ) with these rates as the parame- 
ter a is varied. The only difference between Eq. (fT9"|) and 
Eq. ([7]) with b = —a is an additional constant factor in 
(p~9|) which actually affects the steady state. The analysis, 
however, is similar as we shall see in the next section. If 
we ignore this additional factor of 2 for the moment, the 
amplitude C{a) given by (15TJ) remains positive as long as 
a < i. Since Df — 1/a is the fractal dimension of clus- 
ters, we conclude that the Brownian coagulation with 
input in three dimensions results in a stationary mass 
distribution if Df > 2. The amplitude vanishes when 
Df = 2 thereby questioning that the mass distribution 
reaches a stationary limit in this case. Physically Df = 2 
is realized when clusters are membranes. 

A range of interesting behaviors arises when clusters 
are polymers. Coagulation often generates polymers. For 
instance, this naturally occurs if monomers have two re- 
active bonds. Schematically we can represent a monomer 
as —M-. Pairs of monomers collide to create dimers with 
two reactive bonds, and generally clusters are polymers 
with two reactive bonds: —M = . . . = A/—. In polymer 
physics, we encounter three typical situations |27| . For 
the ideal chain, equivalently a polymer in the so-called 
9 solvent, the polymer is essentially a random walk, so 
Df = 2, or a = h. For a polymer in a good solvent, 
a « |. Based on the discussion above, we expect that 
the mass distribution in this case does not reach a steady 
state [28J . The same is true when polymers are very stiff, 
so that the kernel is given by (fT9")l with a = 1. 

In the following sections we address the challenge of 
understanding driven Brownian coagulation, as charac- 
terized by the Brownian kernel (fT!?)) , in the range of the 
parameter space < a < 1 indicated by the solid line in 
Fig. [TJ The previous discussion suggests that the cases 
a — i a = |, and a = 1 are physically most interesting. 
We first study the stationary case, < a < \ in Sec. Mil 



We then study the non-stationary case, ^ 
Sec. IIVI We leave the marginal case, a = \ 
Sec. [V] since it is the most challenging. 



< a < 1, in 

, to the final 



C. Remarks on numerical simulations 

All numerical simulations presented in this paper were 
done using a variant of the coarse-graining method de- 
scribed in Refs. [29j,|30|. Numerically, one must necessar- 
ily truncate Eq. (|TJ) at some large cluster size which we 
denote throughout by M. The choice of truncation is not 
unique. We chose to truncate by requiring that all clus- 
ters having mass greater than M are removed from the 
system. For all dynamical simulations, care was taken 



Analytic 
Numerics M=10 
Numerics M=10 6 




FIG. 2. Numerical measurements of the amplitude, C, of the 
stationary state as a function of a for the kernel ((TJ) with 
b = —a using two different values for the numerical cut-off, 
M . The theoretical prediction, given by Eq. (|16[) . corresponds 

to M = oo. 



to ensure that simulations were halted before the cluster 
size distribution started to feel the presence of this cut- 
off. We validated the code by comparing the amplitude 
of the stationary state against the exact result, Eq. (|16p . 
The results are shown in Fig. [5] for the kernel (O with 
b = —a. We can see that, our numerical measurements 
are in exact agreement for small values of a. As a ap- 
proaches i, however, there is a considerable finite size 
effect which diminishes as the cut-off, M, is increased. 



III. STATIONARY CASE: < a < \ 



A. Constant kernel 

We begin with the simplest case when the merging rate 
is constant, a = 0. In this setting one can determine both 
the steady state mass distribution and several temporal 
characteristics. (This has been done e.g. [l(J, LLLJ and 
reviewed in 0, Q.) Hence this model sheds the light 
on the emergence of the steady state. Without loss of 



generality, we set Kij = 2 in Eqs. (|2"0|) . This choice 
can be achieved by appropriate rescaling of the density 
and time. The rate equations for the aggregation process 
that proceeds with a constant merging rate and driven 
by input of monomers then read 



dc k 
dt 



E 

z+j— k 



(-',(. j 



2c fc M + 4,i- 



(20) 



The total density of clusters, M = J2k>i c ki is J us t the 
zeroth moment of the cluster size distribution. Summing 




FIG. 3. Numerical calculation of the time evolution of the 
cluster size distribution for the constant kernel. Inset shows 
the data collapse obtained using the scaling given by Eq. (|A2jl 
with a = 3/2. 



Note that c& ~ k 3 / 2 for large k, so that the mass density 
y^ fccfc diverges. This is not surprising since 



Mi =5^fc Cfc (t)=t 

fc>i 



(26) 



due to mass conservation. 

When t is large, the mass distribution Cfc(t) is very 
close to stationary for sufficiently small masses fc « t„ 
while for k 3> fc* the mass distribution is essentially zero. 
The crossover mass fc* is found from 

oo k* k* 

t = ^k fc (t)«^k t ~^r 1 / 2 ~^ (27) 
*;=i fc=i fc=i 

implying that fc* ~ £ 2 . The scaling analysis and numeri- 
cal measurements presented in Appendix [X] suggest that 
fc* ~ t 2 is true for all values of a in the range < a < 1. 



B. Brownian coagulation with input 



Eqs. (l20l) we find that M (t) satisfies 
dM 



dt 



-M' 



1 



so for an initially empty system 

M (t) = tanh(i) 



(21) 



(22) 



Plugging (|2"2")) into the first rate equation ([211)) we obtain 
c±(t) = — 2 tanh(i) ci(£) + 1, from which the density of 
monomers is 



ci(t) = 



1 



f 



cosri (i) 



tanh(i) 



(23) 



Generally, one can solve Eqs. (f2"U|) recursively to give 

Cfc W = ^27T7 / ^'cosh 2 (i') Y, Oi(t')cj(t') (24) 
cosh (t) Jo ^ k 

This formal solution quickly gets very unwieldy. Equa- 
tions (|2"2")) - ([2"3l show that the system reaches a non-trivial 
steady state [31(, Ck{t) — > cv See Fig. [3l 

The determination of Ck from the time-dependent solu- 
tion, (|24p is impractical. A simpler procedure is based on 
the generating function technique introduced in Sec. Ill Al 
We multiply Eq. ([20|) by z k , sum over k and set the left- 
hand side to zero to obtain a stationary state. This trans- 
forms the infinite system (1201) into a quadratic equation 
for the generating function Cq(z): 

= C 2 , - 2C n + z. 



Solving this equation we get Qo( z ) = 1 — yl — z. Ex- 
panding Co (2) we arrive at 



Ck 



1 r(fc-i) 



4tt r(fc + i) 



(25) 
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FIG. 4. Scaling functions (obtained using the scaling given by 
Eq. (|A2| with a — 3/2) multiplied by x 3/2 . Scaling functions 
for a = and a = 1/3 are plotted. The "bump" phenomenon 
at large cluster sizes is clearly evident. 



Consider steady states in a class of models with the 
generalized Brownian kernel (|19|) . Adopting the nota- 
tion of Eqs. (ITU1) and (ITTj) . multiplying (@| by z k and 
summing over &;, we arrive at the following equation for 
the generating functions 

= e a (z)e_ Q (z) + e (z) 2 - e a (z)M_ a - e_ a (*)M a 

-2M e (z) + z. (28) 

Expanding C Q (z), 60(2) and C_ a (z) in the z j" 1 limit 

e„(*) = M a + CT(1 - r - a) (1 - z)^^ 1 + . . . 

C (z) - M + CT(1 - r) (1 - z)^ 1 + . . . (29) 

C_ a (z) =M_ a + CT(1 - r - a) (1 - z^^ 1 + . . . 



and substituting ((29J) into Eq. ([28]) we obtain M a M_ a + 
Mq = 1 in the zeroth order [32J. Matching the leading 
terms in (1 — z) in this expansion, we determine the decay 
exponent to be r = 3/2 and obtain a relation for the 
amplitude C: 



1. (30) 
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Dm this we obtain 
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1 i il 


[(i 


-4a 


2 ) cos(7ra) 



-1/2 



(31) 
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For the classical Brownian coagulation (a = 1/3) we ar- 
rive at the announced prediction © for the tail. Note 
that if one ignores the constant 2 in ([3]), then the fac- 
tor ^/5/23 in © would be replaced by t/5/18 [H[. The 
scaling functions (see Eq. (|A2[1 ) for the constant kernel 
case, a — 0, and the Brownian coagulation case, a = i, 

are plotted in Fig. [4] multiplied by a factor of x? . The 
plateaux for small masses indicate that the scaling ex- 
ponent t = S is very well supported numerically. Fur- 
thermore the predicted values for the amplitudes of these 
plateaux are in reasonable agreement. The discrepancy 
for the case a = \ can be traced to a finite size effect. 
It is interesting to observe that the large mass structure 
of the size distribution is not trivial with a characteristic 
"bump" being clearly visible. We remark that a some- 
what similar bump is observed hydrodynamic turbulence 
as the energy cascade enters the dissipation range. In 
that context, this bump is known as the "bottleneck ef- 
fect" |34| and can be attributed to the depletion of non- 
linear interactions due to the decay of the energy spec- 
trum as it enters the dissipation range. It is tempting to 
suggest, by analogy, that the bump in Fig. [4] is the result 
of the depletion of the reaction rate for large masses due 
to the absence of potential coagulation partners ahead of 
the front. Further investigation would be required to test 
this suggestion quantitatively. In any case, this bump is 
typical for all values of a and plays a very important 
role in the following section where we analyze the case of 
a > i 



IV. NONSTATIONARY CASE ± < a < 1 

A. Stiff Polymers (a = 1) 

The Brownian coagulation of stiff polymers is de- 
scribed by reaction rates 



i 7 
Kij = - + - + 2 

J % 



(32) 



Let us first explicitly demonstrate that the stationary 
state is never reached. For the kernel (|3"2"1) . the first few 



FIG. 5. Time evolution and data collapse for the case of a = 1. 
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FIG. 6. Evolution of the minimum point for the case of a = 1. 



governing equations (UJ) read 

^jr = l-ci(t + M-i + 2M ) 
^ = 2c? - c 2 (|t + 2M_! + 2M ) 
-± = § ClC2 _ c 3 (It + 3M_! + 2M ) 
^± = M ClC3 + 2c 2 - C4 (^ + 4M_! - 



(33) 



2Afn 



where we have taken into account mass conservation (|26[) 
and again adopted the notation (fTTj) to denote moments 
of the cluster size distribution. 

In the long time limit t = Mi 3> Mg > M_i and there- 
fore from (|33p we find the leading asymptotic behavior 
of the densities of light clusters: 

1280 



c\ 



(>2 



t 3 ' 



C3 



54 



C4 



t< 



These partial results reveal the general pattern, viz. 



Bi 



Ck = 



t 2k- 



(34) 



Thus for any fixed mass, the corresponding density ap- 
proaches to zero: lim^oo c& (t) — 0. The formulas ([34|) 
are applicable before the minimum is reached. 

Plugging (|34j) into the governing equations one finds 
that the amplitudes Bk are determined by recurrence 
B x = 1, 



k 

Bk = O Yl K hJ B i B J when k ^ 2 ( 35 ) 

i+j=k 



The sequence B k = 1,4,54,1280,44500, etc. does not 
appear in [35j] and apparently does not admit a compact 
expression that depends only on k. The asymptotic be- 
havior of the amplitudes is simple. We need only one 
relation, BkjBk~\ — k 2 when k 3> 1, which immediately 
follows from (|3"5j) and (|3"2")l . Using this together with (j3~4")l 
one finds that for a fixed large t, the mass distribution 
Ck quickly decreases with k when k < t, a minimum 
Cmin oc e~ 2t is reached at k sa £, and then C/t starts to 
increase. Fig. [5] shows the time evolution of the cluster 
distribution for stiff polymers obtained from numerical 
simulations. The emergence of a minimum is clearly ob- 
served. The predicted linear increase in the position of 
the minimum is well supported by numerics (main panel 
of Fig. [6J. The predicted exponential decrease of c m ; n (t) 
with time is also well supported by numerics (inset of 
Fig. [To]) although the rate of decrease seems to be less 
than 2. 

The behavior of the mass distribution in the k > t 
region is harder to understand than the behavior (j3~4")l in 
the 'boundary layer' region k < t. Before analyzing the 
mass distribution in the k > t region let us look at some 
higher order moments which encode a lot of information 
about the mass distribution. Generally, the moments 
evolve according to rate equations 



dM, 



dt 

dM 2 
dt 



- = ^I>. 



jCiCn 



J,J>1 



= 1 + J^ Kijijacj 



i,J>l 



1 + 3 y, K iJ 2 3 



(36) 



C) ( j 



dt ' ' - ^ "•'-'' J l ° 



dM 4 
dt 



= 1 + 4 ^ K h ji 3 jciCj + 3 ^ Kiji 2 j 2 CiCj 



etc. Specializing (J36I) to the reaction rate (|32|) gives 
dM 



dt 
dM 2 

dt 

dM 3 

dt 

dM± 

dt 



1 - {M 2 + tM„i) 

l + 2(M 2 M Q + t 2 ) 

1 + 3(M 3 M + 3tM 2 ) 

1 + 4M 4 M + 10Af 2 2 + 14iM 3 



(37) 




FIG. 7. Scaling of the moments with time for case a — 1. 

We have analyzed equations (|3"T|) using some guess 
work guided by simulations. The emerging behavior is 
quite simple. The asymptotic of the moment M_i, and 
more generally any moment M p {t) with p < 0, is dom- 
inated by monomers. In particular, the moment M_i 
which appears in Eqs. (|37l) . is given by 



M_! 



(38) 



in the leading order (see red crosses in Fig. [7]) . The mo- 
ments Mp{t) withp > exhibit a simple scaling behavior: 



M p = Apt 2 **- 1 



(39) 



We have Ai = 1 since Mi = t. To find other amplitudes 
we insert (|39p into (|3"T)l and obtain 



3A 2 = 2(A 2 A 
5A 3 = 3(A 3 A 
7A 4 = 4A 4 A ^ 
etc., from which 



-1) 
-3A 2 ) 

10A^ + 14A 3 



A 2 = - 



A' 



A = 3 - 2A 



A, 



A, 



:-!(. 



i,J>l 



i,i>l 



A(1 + 3A) 

40 + 624A 
A 2 (1 + 3A)(1 + 2A) 



(40) 



Generally all A p with p > 2 can be expressed through 
A which therefore remains the sole unknown amplitude. 
Numerics (see inset of Fig. [7j indicate that Ao ~ |. 
These relationships between moments are verified numer- 
ically for the first few moments (the main panel of Fig. [7]). 
The asymptotic behaviors (|38[) -(|39 |) tell us that the 
mass distribution exhibits a scaling behavior outside the 
boundary layer region. If m is the typical mass scale and 
c is the typical density scale, the moments will exhibit 
the scaling behavior 



M p {t) ~ cm p+1 



The consistency between (|39|) and (J4TJ) sets 



r 



m 



• r 



(41) 



(42) 



The same mass scale characterizes the Brownian coagu- 
lation in the situations when the mass distribution ap- 
proaches a stationary limit, namely for the reaction rate 
(fj"9|) with \a\ < 1/2. In this case, Ck ~ A; -3 / 2 when 
k < fe*, and the same argument as in ([27|) gives £;* ~ t 2 . 
The time evolution of the cluster size distribution is 
shown in Fig. [5] and clearly illustrates the predicted sep- 
aration of the size distribution into two populations of 
clusters: a boundary layer near the monomer scale and a 
broader distribution of large clusters. The inset of Fig. [5] 
demonstrates that the large mass part of the cluster dis- 
tribution collapses very well when rescaled using the the- 
oretical scaling given by Eq. (|42|). 

In addition to the scaling cm p+1 contribution to the 
moment M p , there is a contribution coming from the 
boundary layer region where the mass distribution fol- 
lows (|34p . Actually the contribution from the boundary 
layer region is dominated by monomers. The scaling con- 
tribution to the moments M p dominates when p > and 
the monomer contribution dominates when p < 0; at the 
marginal case of p = both contributions are compara- 
ble. These arguments explain (1381) . and more generally 
M p = t- 1 when p < 0, and ([3"9"]l for v > 0. 

Using the scaling form of the mass distribution 

c k (t)=t- 3 F(x), x=^ (43) 

we can express the amplitudes A p via the integer mo- 
ments of the scales mass distribution 



dxF{x) =A -1 



dxxF(x) = 1 



/ dxx p F(x) = A p , p>2 
Jo 



(44a) 
(44b) 
(44c) 



Let us now try to determine the scaled mass distri- 
bution F(x). When i + j = k, the reaction rate (j32|) 
becomes Kij — k/i + k/j and therefore the gain term in 
the Smoluchowski equation (j4|) can be re-written as 



1 fc-i 

— y KijCiCj — K 2_^ Ck-^ 
i+j=k j=l 



J 



(45) 



while the loss term simplifies to 



c k*^2 K k, 3 Cj = c fe 



kM. 



2M ( 



(46) 



The sum on the right-hand side of Eq. (|45[) contains 
the summation over the boundary layer and scaling re- 
gion. The former is dominated by the contribution from 
monomers, while in the scaling region the summation can 
be replaced by integration. Therefore 



fe-i 
k J^ c k - 

3 = 1 



j — ~ fcCiCfc-l + 

J 



fe r dy , x , 

4 i F{y)F(x 



y) (47) 



Inserting Eqs. ([i5 ]) -(|37 jl into (gj) gives 

/* 2A 

t— k) - Ck [k + — 
k_ 



Ck 



(Cfc-l - 

x dy 



o V 



F(y)F(x - y) 



where we have also used M_i = c\ = 1/i, Mq = Ao/t. 
We now replace the difference Ck-i — Ck by the deriva- 
tive and use the scaling form (j43| to find the governing 
equation for the scaled mass distribution 



F' 



1- Xx 



■F 



--F(y)F(x 

V 



y) = (48) 



where, as in Eq. (|40|) . we used notation A = 3 — 2Ao We 
must solve (|48|) subject to 



POO 

I dxxF(x) = l, F(0) = F( 
Jo 



oo) = (49) 



It seems impossible to find an analytic solution of the 
boundary- value problem (J4*5 |) - (|4"9")) . However, we can de- 
termine the asymptotic behaviors and the qualitative 
shape of the scaled mass distribution. For instance, in 
the small x limit keeping the dominant terms simplifies 
to F' = x~ 2 F, from which 



F 



-l/x 



(50) 



In the opposite x — ¥ oo limit, the scaled mass distribution 
decays exponentially 

F~x~ x e-~< x (51) 

Plugging ([ST]) into (l48l) and taking the x — > oo limit we 
arrive at an equation that determines the amplitude 7: 



7 



r^F(y)^ 
Jo V 



Overall, the scaled mass distribution has a bell-shaped 
curve which vanishes exponentially fast in the x — > and 
x — > 00 limits. 

Note that without input the scaled mass distribution 
also has a bell-shaped curve and similar asymptotic be- 
haviors [36[. though the mass scale is different, m ~ t, 
and the scaling ansatz is given by Cfc(i) = t~ 2 F(k/t) in- 
stead of (1431) . The most important distinction from the 
driven case is the lack of the boundary layer region. 



Polymers in Good Solvents: ~ < a < 1 




The same arguments as in Sec. IIV Al show that the 
monomers provide the dominant contribution to the total 
cluster density Mo(t). Thus asymptotically 



M = ci ~t- {2a - l \ 



(56) 



These scaling laws are verified numerically in Fig. [5] We 
remark that there is again a very strong finite size ef- 
fect as a approaches the marginal value of \ so that for 
the physically interesting case of a = | we were unable 
to integrate for times long enough to observe the final 
asymptotic scaling predicted by (|56|) although the trend 
is clearly evident from the behavior of the larger values 
of a. 



MARGINAL CASE OF IDEAL POLYMER 
CHAINS: a = 1/2 



FIG. 8. Verification of scaling of the monomers density pre- 
dicted by Eq. lf56) for values of a in the range i < a < 1. 



The physically interesting situation of polymers in 
good solvents corresponds to a = 3/5, but the follow- 
ing analysis applies to the entire region | < a < 1. The 
procedure is essentially the same as for the a = 1 case. 
The densities of monomers and dimers satisfy 



^=l- Cl (M a + M_ a + 2M Q ) 
at 

^ = 2c 2 - c 2 (2- a M a + 2 a M_ a + 2M ) 



(52) 



Therefore asymptotically we have C\ = 1/M a for 
monomers, c 2 = 2 l+a c\/M a = 2 1+a /M^ for dimers, and 
generally 



Bi 



Ck 



M™- 1 



(53) 



where the amplitudes B^ are determined by recurrence 
B\ = 1 and (|3"B"|) . The mass distribution is described by 
Eqs. (|53p in the boundary layer region. 
The first few integer moments satisfy 

' M/ °=l-(Af 2 + M a M_ a ) 



dt 

dM 2 

dt 

dM 3 

dt 



= l + 2(M 1+Q M 1 _ a + i 2 ) (54) 

= 1 + 3(M 2+a Mi_ a + M 2 _ a M 1+a + 2tM 2 ) 



From these equations we see again the validity of the 
scaling behavior (f39|) which implies Eqs. (J42" |) -(j43" l) . Using 
(|55|) withp = awe conclude that A ~ i 2a_1 and therefore 



C\ ~ t" 



(2a-l) 



(55) 
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FIG. 9. Stationarity of the monomer density for a — 1/2. 



When a = 1/2, the amplitude C(a) given by Eq. ([3T]) 
vanishes indicating that our analysis does not apply to 
this situation. The approach that has led to Eq. (|3"Tj) is 
based on two assumptions: The emergence of a stationary 
solution and an algebraic behavior (Tl"3"|) of the stationary 
solution. Obviously, at least one of these assumptions is 
incorrect when a = 1/2. In this marginal case the ma- 
jor assumption regarding stationarity appears to work, 
while the secondary assumption of the power-law behav- 
ior (fl3|) no longer holds. Since the power-law asymptotic 
is valid for all a < 1/2, one anticipates that the devia- 
tion is merely logarithmic. This suggests the conjectural 
asymptotic 



Ck 



Ck 



-■a/i 



(lnfc)- 



when k — > oo 



(57) 



There are some parallels with a problem in wave turbu- 
lence which we studied previously [37] . In the turbulence 
setting, the corresponding expression for the amplitude 
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diverges rather than vanishes for a particular marginal 
kernel, and we did find that this divergence led to a log- 
arithmic correction. The current problem seems to be 
considerably more difficult as our final asymptotic (|62[) 
involves a product over repeated logarithms. 

Before proceeding we note that the assumption that 
the cluster densities become stationary in the t — > oo 
limit tells us that the moment 



Mi 



fe>i 



('k 



(58) 



must remain finite. Indeed, the rate equation for the 
monomers 



~~dt 



l-d(Mi +M_± +2M ) 



tells us that the moments Mi , M_ 1 and Mn must re- 

2 2 

main finite (of course, the finitcness of the moment Mi 

K 2 

suffices) and the stationary value for the monomer den- 
sity is given by 



1 



Cl 



Mi + M_ i + 2M 



(59) 



The conjectural behavior (|5"T|) of the tail is compatible 
with the finiteness of the moment Mi when ti > 1. Since 

2 

we do not know a-priori that the cluster size distribu- 
tion becomes stationary, we check the finiteness of the 
moments for large times numerically (see inset of Fig. [9]) 
and find fairly strong evidence that the relevant moments 
tend to finite limits as t — > oo in support of the assump- 
tion of stationarity. The stationary value of the monomer 
density is well described by Eq. (|59[) as shown in the main 
panel of Fig. [§] 

To determine the exponent n we use again the gen- 
erating function approach and arrive at ([28]) . The con- 
jectural behavior ([ST)) leads to certain behaviors of the 
generating functions in the z f 1 limit, yet the results are 
inconsistent for all t\ ^ 1. For t\ = 1, the results are 
almost consistent, namely after canceling the dominant 
terms one eventually arrives at a relation of the form 
1 ~ ln(ln(l/2:)). Formally, this is of course inconsistent, 
although the error is much smaller than for any T\ ^ 1. 
Another drawback is that for n = 1 the moment Mi 

2 

diverges, although in a very mild logarithmic manner. 

Thus we should take t\ = 1, although the ansatz (|5"T)) 
with this T\ is still slightly incorrect. A similar prediction 
is obtained by entirely different means in Appendix [B] 
After a bit of trial and error one arrives at an improved 
ansatz 



c fc -> Ck- 3/2 (\nk)- l [\n(\nk)Y 



(60) 



The finiteness of the moment Mi now implies r 2 > 1. 

2 

The singularity analysis of the generating function equa- 
tion (|2"5)l shows that results are inconsistent for all r 2 7^ 1 , 
while for r 2 = 1 the inconsistency is the weakest, namely 
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FIG. 10. (A) Cfe(i) for the case of ideal polymers (a = 1/2), 
(B) multiplied by fc 3/2 , (C) multiplied by fc 3/2 ln(fc) . 



one gets 1 ~ ln3(l/z). (We use notation ln p (a:) for the 
repeated logarithm: By definition lni(x) = ln(x) and 
ln p+1 (a;) = ln(ln p (x)).) Thus we should take r 2 = 1, 



although the ansatz (|60|) with this r 2 is still slightly in- 
correct. Trying 



c fc ^ Cfc~ 3/2 (lnfe)- 1 (ln 2 fcj-^lnafc)- 



(61) 



11 



one again finds that t$ = 1 and the ansatz (|6T|) should 
be correct by another multiplicative factor (ln 4 fc) -1 , and 
then by (Ins k)~ l , etc. In other words 



c k ^Ck-V 2 ]]_(\n p ky 
p>i 



(62) 



A similar subtle behavior characterized by an infinite 
product of repeated logarithms arises in reversible poly- 
merization 38]. In practice, the repeated logarithms are 
essentially undetectable, so at best one can hope to see 



c fc ~fc- 3/2 (lnfc) 



(63) 



These predictions are compared with numerics in Fig. 1101 
Figure [TUtA ) shows the time evolution of the cluster size 
distribution for the marginal case a = \ . The same data 

are shown in Fig. fTOl B) multiplied by a factor of k^ . The 
fact that the tail of the size distribution now looks almost 
flat on this plot indicates that the exponent t = § is cor- 
rect although there is still some weak variation with k in 
the tail which may indicate a logarithmic correction. The 
inset of Fig. rTor B) shows the same data collapsed accord- 
ing to the scaling (|4"2"j) . It is clear that, similar to what we 
found for the non-stationary case in Sec. IIV A| the lead- 
ing bump of the size distribution is very well described by 
scaling whereas the tail is not. Figure [TU^C) shows the 
same data multiplied by an additional factor of log(m) as 
suggested by the theoretical arguments outlined above. 
We concede that it is open to interpretation whether the 
result demonstrates a more convincing plateau for small 
cluster sizes than the corresponding data without the log- 
arithmic correction plotted in Fig. rTOT B). It seems that 
the honest conclusion to draw is that the numerics are un- 
able to give definitive support to the conjectured behavior 
(f6"3")) although the fact that we obtained the same answer 
by completely different means in Appendix [B] gives us 
some degree of confidence that we have obtained the cor- 
rect asymptotic behavior. 



VI. CONCLUSIONS AND DISCUSSION 



polymer aggregates. We find that this behavior falls into 
two classes, depending on the value of a. 

If the exponent of the kernel is in the range < a < | 
then the size distribution for any fixed mass becomes sta- 
tionary for large times (although the typical mass contin- 
ues to grow as t 2 for all times). The resulting stationary 
state is a non-equilibrium stationary state characterized 
by a balance between the generation of new clusters by 
the injection and subsequent aggregation of monomers 
and the loss of clusters via aggregation. The stationary 
state carries a flux of mass through the space of clus- 
ter sizes from small clusters to large. We determined 
that the scaling exponent of this stationary state is al- 
ways | and obtained an exact expression for the ampli- 
tude of this stationary state as a function of a. During 
the course of this analysis, we analytically determined 
the stationary cluster size distribution in the presence of 
source of monomers for the problem of Brownian coagu- 
lation of spherical droplets (a = ,) originally introduced 
by Smoluchowski at the foundation of this field almost 
one century ago. 

If the exponent a is in the range ^ < a < 1, which 
includes the physically relevant cases of stiff polymers 
a = 1 and polymers in an ideal solvent a « | we find 
that the cluster size distribution never reaches a station- 
ary state. Rather the population of clusters splits into 
two populations resulting in a bimodal cluster size dis- 
tribution. In this regime, the size distribution exhibits a 
narrow "boundary layer" of clusters near the monomer 
scale. This boundary layer is separated by a gap, in which 
the cluster size distribution goes exponentially to zero in 
time, from a population of large clusters which continue 
to grow for all times by absorbing small clusters. 

The marginal case, a—\ corresponding to ideal poly- 
mers, turned out to be very difficult to understand theo- 
retically since it exhibits aspects of both behaviors. Our 
analysis suggests that the size distribution in this case 
does become stationary, but the fc -3 / 2 tail acquires a 
logarithmic correction; more work is required to make a 
definitive statement about the final asymptotic state for 
this marginal case. 



To conclude we have presented an extensive analysis of 
the kinetics of Brownian coagulation of polymers in the 
presence of a source of monomers at the mean field level. 
Our study focused on the determining the structure of 
the cluster size distribution for large times as a function 
of the exponent a, the inverse fractal dimension of the 
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FIG. 11. Typical scale, s(t), multiplied by t 2 for several values 
of a. 



Appendix A: Scaling analysis 

The assumption of scaling puts several strong con- 
straints on the dynamics of the cluster size distribution 
which wc summarize here. Consider the continuous ana- 
logue of Eq. ([I]) obtained by replacing, the discrete cluster 
size index, k, with a continuous cluster "mass" , m: 



Cm(*) = ^ 



oo />oo 



dmidm 2 K(mi,m 2 )c mi (t)c rrl2 (t) 



x [S(m — mi — m 2 ) — 5{m — m\) — 6(m — rri2)] 
+ 5(m-l). (Al) 

Let us assume that the solution exhibits scaling. That 
is to say, there exists a monotonically increasing charac- 
teristic cluster size, s(t) and a scaling function of a sin- 
gle variable, F(x), such that the cluster size distribution 
tends to the the scaling form 



c m (t) ^ s(t) a F(x), x 



s(t) 



(A2) 



Throughout this paper, we took the characteristic size to 
be given by the ratio of moments s(t) = M2(t)/Mi(t). 
More precisely, the scaling behavior (JA2I) is supposed to 
emerge in the scaling limit s(t) — > oo and m — > oo with 
the scaling variable x kept fixed. Let us now determine 
the exponent a. Substituting Eq. (|A2|) into Eq. (|Alj) 
we obtain an integro-differential equation for the scaling 
function 

dF [°° 

aF — x—— = —Fix) \ dxi K(x,xi)F(xi) 
dx J 

+ - / / dxidx 2 K(xi 1 x 2 )F(x 1 )F(x 2 )6(x - X\ - x 2 ) 
2 Jo Jo 

and a simple differential equation for the characteristic 
cluster size, 



,A+rt + 2 



(A3) 
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where A = a+b is the degree of homogeneity of the kernel. 
If we further assume that the scaling form contains the 
total mass: 



dt 



mc m (t) dm 



d 
dt 



a a+2 



x F(x) dx 



then we obtain a second equation for s(t): 

s - a -i 
s = 



Fi (a + 2) 



(A4) 



where F% = L x F(x) dx is the first moment of the scal- 
ing function. Comparing Eqs. (|A3|) and (1A4[) we obtain 
the exponent a: 



A + 3 



(A5) 



Equations (JA3I) and (IA5|) require that the typical size 
grows as a power law for large time: 



s(t) ~ £T^A 



as t 



oo. 



(A6) 



For the case of Brownian coagulation, A — 0, the as- 
sumption of scaling together with the assumption that 
the mass is concentrated in the scaling part of the size 
distribution, requires that s(t) ~ t 2 . This prediction is 
in good agreement with numerics. See Fig. [TT] for some 
representative numerical results for the generalized Brow- 
nian kernel, Eq. (|19|l . for several values of a. 

Furthermore, if the size distribution becomes station- 
ary (and non-zero) as t — > oo, then Eq. (|A2|) and Eq. (|A5[) 
require that the scaling function must be algebraic for 
small values of x in order to cancel the time dependence: 



F(x) 



as x 



0. 



(A7) 



Again for the case of Brownian coagulation, assuming 
scaling leads us to expect that if we have a stationary 
state it must scale as k~ 3 ^ 2 for small cluster sizes. 



Appendix B: An alternative derivation of the 
stationary state 

An alternative way of obtaining the stationary state 
amplitudes given by Eq. (IT^|) or Eq. ([3~Tj) was outlined in 
yjj. We summarize this method here since it provides 
some insight into what happens in the marginal case, 
a = i. If we assume that the size distribution behaves 
algebraically in the large m limit, c m = C ' mT T , then 
after some re- arrangement, (|A1|) can be written as 



C 2 



OC /"OO 



dm i dm 2 



(Bl) 



\K(mi, 7712) (mi 7712) r 6(m — mi — 7712) 
—K(m, mi) (mmi)~ T 5(m 2 — m — mi) 
— K(m, m-i) (mm2)~ T 5(mi — m 2 — m)] 



The right-hand side contains three integrals and we 
now apply the following changes of variables 



(7711,7712) 



(mi,m 2 ) 



777777, 777 



777J, 777,- 



777 mm-. 



777-1 777-, 



to the second and third integrals, respectively. After per- 
forming some algebra and taking advantage of the fact 
that the kernel is a homogeneous function of its argu- 
ments, we obtain 



1 



- m 



a+i-2t r a 



C 2 I(t) 



(B2) 



where 



1 /.1 



I(r) =11 
Jo Jo 



dfiidfx 2 K(^i,fi 2 ) (p-ifJ-2)' 



(i-^r x - 2 -^- x - 2 )s(i-vi-^). 

It is clear that a stationary state is obtained if the expo- 
nent r is given by r = (A + 3)/2. The conservation law 
expressing the conservation of mass, 



d(mc m ) d,J„ 



dt 



dm 



(B3) 



defines a flux of mass, J m , through mass scale 777. Mul- 
tiplying Eq. (|B2j) by 777 and integrating therefore allows 
us to express this mass flux as a function of m for any 
value of the exponent r: 

7 / x lC 2 m x + 3 - 2T I(T) 

J "' (T) = -2 A + 3-2r ■ (B4) 

From Eq. (|B3I) . we see that the flux should become in- 
dependent of 777 (and equal to 1 in our units) in order 
to have a stationary state. Thus, we can evaluate the 
amplitude C as 



C 



lim 



/2(2r-A-3) 



h+l V 777 A + 3 " 2t /(t) 

1 dl(r) V 5 



dr 



(B5) 



the latter step resulting from the use of l'Hopital's rule 
to evaluate the initial undetermined expression. For the 
general Brownian kernel, Eq. (flU|) . we obtain 



C(a) = 



1 



d/i 



(B6) 



2j [MI-m)] 372 

[v a (i-vT a + (i-ri>- a + 2] 

[Mln(/i) + (l- M )ln(l- M )] 



We have not succeeded in analytical computation of the 
integral in Eq. (JB6I) and demonstration that the result is 
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FIG. 12. Comparison of expressions (|B6|) and (|31|l for the 
amplitude of the stationary state. 



identical to (|3Tj) which was obtained via the generating 
function route. Numerical integration (Fig. rT12j) indicates 
that the two expressions are indeed identical. 

While this approach has yielded an answer which we 
already established before, it has the advantage of work- 
ing for any homogeneous kernel. Furthermore, it is very 
helpful in providing physical insight about what happens 
for the marginal case, a = 5, when C(a) vanishes. If we 
formally repeat the preceding calculation taking into ac- 
count the presence of a small mass cut-off at the monomer 
scale, we would obtain [39| 



C(a, m) = 



df 1 



(B7) 



[m(i - M)] 3/2 

[/i]nGu) + (l-/i)]n(l-AO] 



Analysis of the divergence of the integrand in Eq. (|B7|) 
at the endpoints of the region of integration indicates 
that this divergence is integrable if a < \ . For such ker- 
nels, therefore, this integral becomes independent of the 
monomer cut-off when we consider m to be much larger 



than the monomer scale and the amplitude of the sta- 
tionary state is given by a universal constant. As we 
have seen above, this amplitude also characterises inde- 
pendence of the flux of mass through any given mass 
scale, m, in the stationary state. This calculation tells 
us that when we consider masses much larger than the 
monomer mass, if a < i there is no contribution to this 
flux from aggregation with monomers. The mass transfer 
is thus local in the mass space in the sense that mass is 
transferred primarily through the aggregation of compa- 
rable sized clusters. This is contrast to nonlocal models 
where mass is transferred primarily through the aggre- 
gation of large masses with small masses of the order 
of the monomer scale. From this perspective, the addi- 
tion model [40l | is the most extreme example of non-local 
transfer. The corresponding analysis of the more general 
class of models given by Eq. ([7]) shows that the mass flux 
is local in the stationary state provided \a — b\ < 1. This 
is the physical origin of the region of regular stationary 
behaviour in Fig. [1] 

When a = |, the integrand in Eq. (|B7|) diverges as 
— 21n(/x)//x at the lower cut-off which integrates up to 
give a logarithmic dependence of Cf^,™) on m as m 
gets large: 



C(^m) 



ln(m) 



(B8) 



This suggests that the stationary distribution for ideal 
polymers acquires a logarithmic correction: 



Ck 



2k-hn(k)- 1 



as suggested by the previous analysis of Sec. [V] Here it 
should be noted that the assumption that C(a, m) be- 
comes independent of m as m becomes large and thereby 
allowing it to be taken outside of the integral in Eq. (|B1I) , 
remains marginally inconsistent as m grows. Attempting 
to correct this inconsistency by incorporating this loga- 
rithmic correction into Eq. (|B1[) seems likely to lead again 
to the repeated logarithms of Sec. [V] For this reason, al- 
though our numerics are not definitive on this matter, we 
believe the evidence is in favour of stationary behaviour 
for the case of ideal polymers but with logarithmic cor- 
rections to the fc -3 / 2 scaling of the size distribution. 



